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VARIABILITY OF THE VELA PULSAR- WIND NEBULA OBSERVED WITH CHANDRA 
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ABSTRACT 

The observations of the pulsar-wind nebula (PWN) around the Vela pulsar with the Advanced CCD 
Imaging Spectrometer aboard the Chandra X-ray Observatory, taken on 2000 April 30 and November 30, 
reveal its complex morphology reminiscent of that of the Crab PWN. Comparison of the two observations 
shows changes up to 30% in the surface brightness of the PWN features. Some of the PWN elements 
show appreciable shifts, up to a few arcseconds (~ 10 16 cm), and/or spectral changes. To elucidate the 
nature of the observed variations, further monitoring of the Vela PWN is needed. 

Subject headings: supernova remnants: individual (Vela) — pulsars: individual (PSR B0833-45) — 
X-rays: individual (Vela pulsar-wind nebula) 



1. INTRODUCTION. 

X-ray and radio observations show that many of young 
pulsars arc enveloped by compact nebulae emitting a non- 
thermal spectrum (e.g., Gaensler 2000). It is commonly 
accepted that this emission is synchrotron radiation from 
a relativistic pulsar wind shocked in the ambient medium. 
The most famous and best studied example is the Crab 
pulsar- wind nebula (PWN). High- resolution observations 
of the Crab PWN have revealed its remarkably complex 
morphology in radio, optical and X-rays (Bitenholz & Kro- 
nberg 1992; Hester et al. 1995; Weisskopf et al. 2000, and 
references therein). Moreover, the multi-year studies have 
shown complex variability of the Crab PWN. For instance, 
the famous "wisps" , discovered in optical observations, 
show changes in their structure and brightness (Scargle 
1969; Hester et al. 1995), while the bright "knots" change 
their location and appearance between images separated 
by only six days (Hester 1998). Recent observations with 
the Chandra X-ray Observatory (CXO) have shown varia- 
tions on a timescale of several weeks (D. Burrows & J. Hes- 
ter, personal communication). Greiveldinger & Aschen- 
bach (1999) analyzed five ROSAT observations of 1991-97 
and found monotonic changes in the surface brightness, 
some increasing and others decreasing, at a rate of ~ 2% 
yr . Studying the PWN variability offers a unique op- 
portunity to understand the structure and dynamics of 
the relativistic pulsar winds, elucidate the mechanisms of 
PWN formation, evolution and interaction with the ambi- 
ent medium, and establish the properties of the relativistic 
plasmas in PWNe. Thanks to the superb angular resolu- 
tion of CXO, a unique opportunity for such investigations 
is provided by the Vela PWN. 

Observations of the Vela pulsar with Einstein (0.1-4 keV 
band) have shown that it is embedded in a "kidney-bean" 
nebula of ~ 2' size, which emits a power-law spectrum 
with a photon index 7 = 1.7 ± 0.2 (Harnden et al. 1985). 
The Vela PWN was further studied in soft X-rays with 
EXOSAT (Ogelman & Zimmcrmann 1989) and ROSAT 
(Ogelman, Finley, & Zimmermann 1993; Markwardt & 
Ogelman 1998), and it has also been detected at higher X- 
ray energies with the Birmingham Spacelab 2 (2.5-25 keV; 



Willmore et al. 1992) and Compton Gamma-ray Observa- 
tory (44-370 keV; de Jager, Harding, & Strickman 1996). 
Ogelman & Koch-Miramond (1989) claimed optical detec- 
tion of a diffuse nebula around the pulsar, with a size of 
40"-90" and V magnitude of 16-17, which has not been 
confirmed by later observations (Mignani et al. 2001). No 
extended radio emission, at a level of 60 //Jy per 20" x 10" 
beam at A = 6 cm, was found at the site of the X-ray 
PWN, although a ridge (3' x 0.'7) of highly polarized radio 
emission was detected at ~ 2' NE of the pulsar (Bitenholz, 
Frail, & Hankins 1991). 

First CXO observations of the Vela pulsar and its PWN 
(Pavlov et al. 2000; Helfand, Gotthelf, & Halpern 2001; 
Pavlov et al. 2001; Kargaltsev et al. 2001) have shown 
a spectacular fine structure which resembles the Crab 
PWN — arcs, jets, knots, and diffuse cometary tails 
(see Fig. 1). The symmetry axis of the PWN image 
(P.A.= 307° ± 2°), which can be interpreted as the pulsar 
spin axis, is nearly parallel to the direction of the pul- 
sar proper motion (P.A.= 305°, Bailes et al. 1989). This 
suggests that the "natal kick" of the pulsar was directed 
along the rotation axis of the neutron star progenitor (e.g., 
Lai, Chernoff, & Cordes 2001). The spectra of PWN ele- 
ments are power laws of different slopes, which vary from 
7 ~ 1.3 for the NW jet and inner arc to 7 ~ 1.7 for the 
outer diffuse nebula (Kargaltsev et al. 2001). The X-ray 
luminosity of the whole nebula in the 0.1-10 keV range 
is L nch = (6.0 ± 0.5) x 10 32 erg s" 1 (for d = 250 pc - 
Cha, Sembach, & Danks 1999), only ~ 10 -4 of the pulsar 
spin-down luminosity E — 6.7 x 10 36 erg s _1 (vs. 0.045 
for the Crab PWN). 

First indication on the Vela PWN variability was re- 
ported by Helfand et al. (2001) who analyzed two images 
obtained with the CXO High Resolution Camera on 2000 
January 20 and February 21. They noticed a 5% brighten- 
ing of the outer arc and suggested that it may be connected 
with the large pulsar glitch of 2000 January 16 (Dodson, 
McCulloch, & Costa 2000). In this Letter we report new 
results on X-ray variability of the Vela PWN obtained with 
the CXO Advanced CCD Imaging Spectrometer (ACIS; 
Garmire et al. 2001). 
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2. OBSERVATIONS 

Two ACIS observations of the Vela pulsar and its PWN 
were carried out on 2000 April 30 and November 30 with 
exposure times 10,577 s and 18,851 s, respectively. In 
both observations the target was imaged on the back- 
illuminated chip S3. To image the whole PWN on one 
chip, the pulsar was offset from the ACIS-S aimpoint by 
— 1.'5 along the chip row. To reduce pile-up, we used 1/2 
subarray, with a frame time of 1.5 s. For the analysis, we 
used the pipe-line processed Level 2 data (versions 13.2 
and 12.1 for the first and second observations, respectively) 
and the CIAO software, v. 2.0. A detailed analysis of the 
PWN spectra and morphology will be presented elsewhere 
(Kargaltsev et al. 2001). Herein we concentrate on com- 
parison of the two observations. 

To compare the PWN structures, it is particularly im- 
portant to co-align the images properly. Since no X-ray 
point sources, other than the piled-up pulsar, are detected 
in the 8' x 4' field of view, we have to rely upon the aspect 
reconstruction which translates the actual event positions 
on the chip into sky coordinates. For the processing soft- 
ware versions used, the rms aspect offset (error of celestial 
location) 1 is smaller than 0'.'6, although in some cases the 
offset can be as large as 2". To estimate the errors in 
our case, we found centroids of radial distributions of pul- 
sar counts, which correspond to a = 8 h 35 m 20!628, 5 = 
-45 o 10'34'/93, and a = 8 h 35 m 20!614, 6 = -45°10'35'.'26, 
for the first and second image, respectively. The difference 
of these positions, 0'.'36, is comparable with our estimated 
centroiding error, ~ 0'.'5 (about one ACIS pixel). 

The images of the central part of the PWN are shown in 
the upper and lower panels of Figure 1 for the first and sec- 
ond observation, respectively. The white contours in the 
upper panel enclose PWN regions chosen for the compar- 
ison with the second observation. The regions numbered 
1 through 7 enclose well-defined bright features such as 
the outer arc (2a+2b) with a brightened spot (1) at its 
apex, NW jet (3), inner arc (4a+4b), SE jet (5), NE knot 
(6), and SW knot (7). In addition, to examine the vari- 
ability in the diffuse emission, we define regions 8a and 
8b. The white contours in the bottom panel are plotted at 
the same positions as those in the upper panel. A visual 
comparison of surface brightnesses within the white con- 
tours immediately shows substantial differences between 
the two images — e.g., the SW jet is brighter, the spot is 
dimmer, and the NE knot is invisible in the second image. 
Moreover, we see that some elements appear at different 
locations — the SW knot moved 2" westward, and the 
outer arc shifted by, on average, 2" north-west from the 
previous positions. Therefore, we defined new positions of 
six PWN elements in the second image with blue contours 
which enclose the same areas as (and are congruent to) the 
corresponding white contours. To quantify the changes in 
the surface brightness, we measured the numbers of counts 
per unit area per unit time within the white and blue con- 
tours for the first and second observations, respectively 
(see Fig. 2), and found that the changes are indeed quite 
significant for some elements — e.g., 32% ± 4% for the SE 
jet, -26% ±2% for the spot, and -19% ±2% for the outer 
arc (the uncertainties are la statistical errors), while they 

1 see http://asc.harvard.edu/mta/ASPECT/cel_loc/cel_loc.html 

2 see http://asc.harvard.edu/clao/wrkshp/esal.pdf 



are comparable with statistical fluctuations for the oth- 
ers (e.g., for the shifted inner arc). To visualize the spa- 
tial distribution of the brightness/morphology change, we 
scaled the images to the same exposure time (10,577 s), 
adaptively smoothed them with CIAO task csmooth, and 
subtracted the first image from the second one. The dif- 
ference image (Fig. 3) clearly shows the displacements of 
the arcs and the SW knot, brightening of the SE jet, and 
dimming the spot and the NE knot. It also shows some 
brightening (dimming) of the diffuse emission at the NE 
(SW) outskirts of the PWN and the large displacement, 
w 8" toward SW, and shape variation of the filamentary 
structure in the upper right corner of the image (75" NW 
from the pulsar). 

To evaluate systematic errors of the brightness changes, 
one should take into account that, although the pulsar 
was imaged at almost the same location on the chip in 
the two observations, the PWN regions are imaged on dif- 
ferent sites because the roll angles were different (260? 3 
and 55?2). Since the CCD quantum efficiency varies over 
the chip 2 , this may lead to artificial changes of brightness. 
To examine this effect, we employed two sets of on-orbit 
calibration exposures (28 ks and 35 ks, at dates close to 
those of our observations), during which the ACIS chips 
were illuminated by on-board radioactive X-ray sources. 
For each of the calibration datasets, we measured the sur- 
face brightnesses of eight domains on the chip where four 
PWN elements (regions 1, 2b, 4a and 5) were imaged in 
the first and second observations of the Vela PWN. The 
differences of the brightnesses of the domains correspond- 
ing to the same PWN elements, both in the separate cali- 
bration datasets and between the two datasets, are about 
2%-3%, comparable to the statistical errors. This puts an 
upper limit of 3% on both spatial and temporal variabil- 
ity of quantum efficiency within the chip area (<~ 1' x 1') 
where the bright PWN elements were imaged. As an in- 
dependent test, which additionally accounts for the effect 
of dither on the exposure times, we constructed exposure 
maps (e.g., Davis 2001) for the two PWN observations, 
for an energy of 1 keV close to the maxima of the count 
rate spectra. Inspection of these maps shows that the ef- 
fective area varies by about 20% over the whole 8' x 4' 
ficld-of-view, but these large variations are associated with 
the boundaries between the four chip nodes. Since the 
PWN elements under investigation are all imaged within 
the same node 1, the variations are much smaller, w 2%. 
Thus, the large, ~ 20%-30%, brightness changes we de- 
tected is not an instrumental effect — they characterize 
real PWN variations. 

We have also measured the mean surface brightnesses in 
the two observations within a 30" radius circle (centered 
on the pulsar) . To take into account the above-mentioned 
nonuniformity, we divided the original images over the ex- 
posure maps and found a change of —3.0% ± 0.5%. Al- 
though this difference looks statistically significant, its 
magnitude is comparable with possible systematic errors 
(e.g., caused by inaccuracy of the monoenergetic exposure 
maps we used). 

Comparison of the two observations allows one, in prin- 
ciple, to check whether the pulsar luminosity has changed 
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in 7 months. A change of the luminosity could be caused 
by the strong glitch of 2000 January 16 because glitch ef- 
fects can manifest themselves on time scales from weeks 
to years, depending on the depth where the energy release 
occurred. We measured the radial count distribution of 
the piled-up pulsar image and found that its very central 
part, within a 1" radius, became brighter by 30% ± 5%, 
whereas the brightness did not change at larger radii. In 
obvious contradiction with this result, the difference in the 
pulsar countrates, —8% ±9%, estimated from the "trailed 
images" (one-dimensional images formed during the frame 
read-outs) is statistically insignificant. Since the trailed 
images do not suffer from the pile-up, we consider the lat- 
ter result more reliable. The apparent variation of the 
central part of the pulsar image is likely caused by a small 
fluctuation of observing conditions (e.g., focal plane tem- 
perature), which may lead to a considerable change in the 
nonlinear regime associated with strong pile-up. It should 
also be mentioned that the pilcd-up images are very asym- 
metric. In particular, in both observations we see a "blob" 
of 0'.'7 radius at a distance of l'/5 from the center of the 
pulsar image (see Fig. 3), at an angle of about 22° from 
the Z + axis of the spacecraft. The blob is most likely due 
to a tilt between the telescope mirrors in the innermost 
shell (Jerius et al. 2001). It can hardly be seen when there 
is no pile-up because it is much fainter than the core of 
the point source image, but it appears quite bright in a 
strongly pilcd-up image because the pile-up reduces the 
core brightness. 

In addition to the shifts and brightness changes, we have 
examined the spectral changes of the PWN elements. To 
characterize the spectra, we use the hardness ratio ft.1.3, 
defined as the ratio of counts above and below E = 1.3 
keV 3 . For the comparative analysis, this empirical quan- 
tity is more convenient than spectral fitting parameters 
because it does not depend on the spectral model, and 
does not suffer from uncertainties of the detector spectral 
response. We see some evidence of hardness changes, up 
to +29% ± 8% in the SE jet (see Fig. 2), but their sta- 
tistical significance, < 3.6a, is not as high as that of the 
brightness changes. 

3. DISCUSSION. 

The comparison of the two ACIS observations taken 7 
months apart shows that the Vela PWN is by no means 
static — we have detected considerable brightness changes 
and/or shifts of some PWN elements, which are likely ac- 
companied by spectral changes. With only two observa- 
tions carried out, we cannot trace these variations to es- 
tablish their time scales and velocities. In particular, it is 
hard to conclude whether the apparent shifts of the arcs 
and the SW knot are associated with steady motion of 
matter during the 7 months or are manifestations of wave 
phenomena, or some instabilities, in the post-shock rela- 
tivistic plasma. The fact that the NE knot has disappeared 
and the other PWN elements have changed their bright- 
ness implies that the changes are more complicated than 
simple motion of plasma inhomegeneities. 

The remarkable similarity of the Vela and Crab PWN 
morphologies suggests that their variabilities are driven by 

3 Being denned as a linear function of pulse-height amplitude (sec http 
E corresponds to the most probable photon energy for a recorded Pt 



similar physical processes. The main distinction between 
the two PWNe is in their sizes and energetics. In par- 
ticular, the physical size of the Vela PWN, <~ 0.1 pc at 
d = 250 pc, is an order of magnitude smaller than that of 
the Crab PWN, in rough correspondence with the scaling, 
r s oc E 1 / 2 , for the shock radius in an ambient medium 
with a given pressure. This correspondence indicates that 
typical pressures, p ~ E/Ancr 2 ~ 10 -9 erg cm -3 , are 
of the same order of magnitude, which is in agreement 
with the close values of their equipartition magnetic fields, 
B ~ 10~ 4 G (Kargaltsev et al. 2001). The similarity of 
physical parameters of the relativistic plasmas in the two 
PWNe suggests that typical velocities of MHD waves, pre- 
sumably responsible for some of the observed variabilities, 
are also similar. In a relativistic magnetized plasma these 
velocities depend on the mean energy density e, magneti- 
zation parameter a = B 2 / (Aire), and the angle 9 between 
the wavevector and the magnetic field (e.g., Gedalin 1993). 
For the conditions expected in the Crab and Vela PWNe, 
they vary from very low values, w (c/2)v / 3o : | cos#|, for the 
Alfven and slow magnetosonic waves at a <C 1, to ~ (0.6- 
0.7)c for the fast magnetosonic wave at a ~ 1. The fastest 
speeds observed for the Crab PWN wisps, about 0.5c (Hes- 
ter 1998), are thus close to the upper end of the expected 
velocity range. If plasma perturbations can propagate 
with similar velocities in the Vela PWN, the observable 
shifts ~ l"-2" can occur in just few days, so that the 7 
months between our observations may not be a represen- 
tative time period to investigate the variations. 

The different behavior of the Vela PWN elements sug- 
gests that various processes are responsible for the PWN 
variability — e.g., the apparent motion of the arcs can 
be caused by wave processes, while the brightness varia- 
tions of the jets can be associated with large-scale inhomo- 
geneities of the polar outflows. A variety of mechanisms 
have been invoked to explain the Crab PWN variability 
- e.g., formation of MHD shock waves in slightly inho- 
mogeneous wind streams (Lou 1998), instabilities driven 
by synchrotron cooling in the flow (Hester 1995), nonlin- 
ear Kelvin-Helmholtz instabilities in the equatorial plane 
of the shocked wind (Begelman 1999), or cyclotron insta- 
bilities of ion rings (Spitkovsky & Arons 2000). To investi- 
gate the actual roles of such processes in the Vela PWN, its 
behavior should be monitored with CXO — other X-ray 
missions do not have sufficient angular resolution, while 
observations outside the X-ray band are extremely diffi- 
cult, if possible at all, because of the PWN faintness in 
the optical and radio. Moreover, the Vela PWN is more 
suitable for studying brightness and, particularly, spectral 
variabilities with the CXO ACIS because it is not as bright 
as Crab, and its ACIS images do not suffer from pile-up 
which strongly complicates the data analysis. Thus, fu- 
ture CXO observations of the Vela PWN will provide the 
unique opportunity to elucidate the physical properties of 
the relativistic plasmas in PWNe. 

We are grateful to Leisa Townsley, George Chartas and 
Slava Zavlin for useful discussions. This work was sup- 
ported by SAO grant GO1-2071X and NASA grant NAS8- 
38252. 

asc .harvard. edu/cal/Links/Acis/acis/Cal_prods/gaincti/06_19_00/cti .html 
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Fig. 1 — ACIS-S images of a central part (57" X 55") of the Vela PWN of 2000 Apr 30 (top) and Nov 30 (bottom). The pixel size is 0'.'492. 
The white contours in the top panel define the PWN elements in the first image. In the bottom panel, the white contours correspond to the 
contours in the top panel, while the blue contours demonstrate the displaced nebular elements in the second observation. The brightest spot 
is the Vela pulsar. 
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Fig. 2. — Surface brightness vs. hardness ratio hi. 3 for different PWN regions for the two ACIS-S observations. For the second observation, 
the values for the regions 1-4 were calculated within the blue contours in the bottom panel of Fig. 1. (The differences would be even larger 
if exactly the same regions were used for the two observations.) 



Fig. 3. — The difference image (142" X 142") of the Vela PWN. The color scale is in counts ks _1 arcsec -2 . The dashed rectangle corresponds 
to the size of images in Fig. 1. The blue and red structures within the 3'.'5 circle at the center correspond to the "blob" in the pulsar image 
in the first and second observation, respectively (see text for details). The blue and yellow linear structures are the pulsar trailed images. 



